基于 DBN 和 多 传感器 数据 分 解 的 核电 机 组 振动 故障 检测 方法 


摘要 : 由 于 仅 提 取 了 核 


数据 i 


行 平滑 计算 与 融合 处 理 ， 在 DBN 的 作 


包机 组 振动 信号 的 单一 特征 
了 基于 DBN 和 多 传感器 数据 分 解 的 核电 机 组 振动 故障 检测 方法 。 
] 下， 提取 出 核 ! 


， 导 致 核 


获取 核 


包机 组 振动 故障 检测 方法 的 检测 效果 不 佳 。 因 此 ， 设 计 
包机 组 振动 信号 数据 ， 对 获取 的 多 传感器 


机 组 振动 信号 的 多 个 特征 ， 计 


不 同 特征 的 敏感 性 


指数 和 模糊 焙 ， 对 振动 信号 特性 进行 分 析 ， 构 建 对 应 的 振动 故障 检测 模型 ， 求 解 核电 机 组 振动 故障 信号 。 实 验 结果 


表明 ， 在 该 方法 的 实际 应 用 中 ，AUC-ROC | 
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| 线 面 积 更 接近 于 1， 检 测 效 果 更 好 。 
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Vibration fault detection method for nuclear power units based on DBN and multi-sensor data decomposition 


Abstract: Due to only extracting a single feature of the vibration signal of the nuclear power unit, the detection effect of the 


vibration fault detection method for nuclear power units is poor. Therefore, a vibration fault detection method for nuclear 


power units based on DBN and multi-sensor data decomposition was designed. Obtain vibration signal data of nuclear 


power units, smooth and fuse the obtained multi-sensor data, extract multiple features of the vibration signal of nuclear 


power units under the action of DBN, calculate the sensitivity index and fuzzy entropy of different features, analyze the 


characteristics of the vibration signal, construct a corresponding vibration fault detection model, and solve the vibration 


fault signal of nuclear power units. The experimental results show that in the practical application of this method, the 


AUC-ROC curve area is closer to 1, and the detection effect is better. 
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0 引言 


核电 机 组 振动 故障 的 检测 和 预防 十 分 


要 。 然 而 ， 


障 检测 的 准确 
提取 和 模型 


传统 的 故障 检测 方法 主要 基于 单一 传感器 进行 数据 采 


集 和 分 析 ， 
主 会 出 现 误 判 或 漏 检 的 情况 。 


~ 


文献 [1 利用 传感器 实时 采集 医 ; 


这 种 方法 在 面 对 复 杂 多 变 的 振动 


言 号 时 ， 往 


并 对 采集 到 的 信号 数据 


行 去 品 ， 再 从 中 提 


故障 相关 的 特征 


， 在 ThinkPHP 框架 的 作用 


类 器 模型 


上 4， 由 此 检测 出 振动 信号 中 的 故障 信 


订 器 械 的 振动 信号 ， 


理 


取出 与 振动 
下 ， 构 建 分 
号 ， 能 够 实 


时 采集 并 处 理 医 疗 器 械 的 振动 数据 ， 及 时 发 现 潜 在 的 故 


障 , 但 其 


需要 实时 处 理 大 量 的 振动 


文献 [2] 对 采集 的 振动 信号; 
绕组 松动 故障 相关 的 特征 ， 
数 ， 对 其 


故障 特征 信息 ， 提 高 了 检测 和 


言 号 , 检测 效率 较 低 。 
行 频谱 分 析 ， 从 
由 此 计算 变压器 的 频 响 
变压器 中 的 故障 进行 分 析 和 讨论 ， 根 据 提 取 的 


中 提取 出 


与 
函 


吉 果 的 可 靠 性 ， 


应 用 中 对 传感器 的 灵敏 度 要 求 较 高 ， 导 致 其 
生 较 差 。 文 献 [3] 将 采集 的 大 型 涡轮 机 械 i 
进行 滤 ; 
自动 运行 包 络 线 的 作用 下 ， 确 定 涡轮 机 械 在 
并 实时 监测 涡轮 机 械 的 瞬 态 数 
据 ， 一旦 出 现 数据 超过 取 值 范围 
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实时 
为 基础 ， 去 除 其 中 的 异常 值 ， 


各 种 参数 的 取 值 范围 ， 


旦 其 在 实际 


再 对 其 


检测 方法 的 
运行 数据 作 
皮 处 理 ， 在 
正常 状态 下 


， 则 认为 当 


故障 ， 通 过 


自动 生成 运行 包 络 线 和 实时 监测 
极 大 地 减少 了 操作 难度 ,但 其 仅 适 用 


于 大 型 


检测 ， 适 用 范围 


较 窄 。 文 献 [和 实时 采集 风力 


前 存在 瞬 态 
屏 态 数据 ， 
涡轮 机 械 的 
发 电机 的 运 


行 数据 ， 并 对 运行 数据 i 


的 支持 下 ， 从 中 提取 与 振动 故障 相关 的 特征 


行 清洗 ， 在 自动 特征 学 习 


法 


参数 ， 由 此 


构建 对 应 的 故障 检测 模型 ， 通 过 调整 模型 参数 ， 对 已 有 


的 故障 样本 进行 学 习 ， 实 现 对 风力 发 电机 振 


动 故障 的 检 


测 ， 能 够 准确 检测 出 风力 发 


包机 的 振动 故障 


， 提 高 了 故 


想 的 结果 ， 
在 以 往 丰 
传 感 数据 分 解 的 核 
器 数据 分 解 的 前 提 下 ， 将 核 ! 
合 ， 创 新 性 地 利用 DBN 对 核电 机 组 振动 信号 的 特性 i 
行 分 析 ， 从 而 构建 对 应 的 振动 故障 检测 模型 ， 提 高 了 核 
电机 组 振动 故障 检测 的 + 
现 潜在 的 故障 隐患 。 
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| oa\eh 


检测 结果 存在 不 稳定 性 。 
究 的 基础 上 ， 本 文 设计 了 基于 DBN 


性 , 但 其 需要 大 量 的 振动 数据 来 进行 特征 
对 于 数据 量 不 足 的 情况 可 能 无 法 获得 


和 多 


E 确 性 和 可 靠 性 


1 核电 机 组 振动 故障 检测 方法 设计 
1.1 基于 多 传感器 数据 分 解 的 核电 机 组 振动 信号 融合 


核电 机 组 振动 信号 能 够 真实 


行 状态 5 


， 对 划 


处 理 。 


与 其 


复杂 性 和 多 变 ， 且 在 采集 过 程 中 


中 可 
他 数据 特定 不 同 ， 核 


型 不 同 ， 导 致 采集 的 振动 信号 也 存在 一 定 差异 四。 
在 对 核电 机 组 振动 故障 进行 检测 


机 组 振动 故障 检测 方法 。 在 多 传 感 
电机 组 振动 信号 数据 进行 融 


， 有 助 于 及 时 发 


地 反映 核电 机 组 的 运 
能 出 现 故障 问题 ， 进行 及 时 维护 和 
包机 组 的 振动 信号 具有 
j 的 传感器 类 


因此 ， 
时 ， 需 要 先 利 用 不 同类 


别 的 传感器 对 核电 机 组 振动 信号 进行 采集 ， 从 而 获取 多 


传感器 数据 门 。 


其 获取 的 多 传感器 数据 如 下 所 示 : 


X=C,+D,cos(wt+p,) (1) 


传感器 数据 ， C; 表示 振动 信号 在 不 同 传感器 的 直 


量 ，D, 表示 不 同 传感器 采 


在 上 述 


公式 中 ，x 表示 获取 的 术 


机 组 振动 信 


号 多 


流 分 


集 到 的 振动 信号 交叉 流 


量 最 


呀 


大 值 ，w 表示 核电 机 组 在 不 同时 间 段 的 角 频率 ，t 表示 


振动 信号 采集 的 总 时 长 ， 有 表示 不 同 传感器 的 相位 角 。 


为 提高 上 述 获 取 的 多 传感器 数据 的 质量 ， 需 要 对 其 进行 
平滑 计算 。 其 具体 计算 过 程 如 下 所 示 : 


B(y)=B| FRcos(wt+p)| 
E -eos(wetB)+ 走 


S(x)=B|E*B(y) Et | 


[a 


在 上 述 公式 中 ，B(y) 表 示 上 述 获取 的 多 传感器 数 


据 的 数学 期 望 值 ， 书 表示 多 传感器 数据 的 波形 参数 ，E 
表示 在 获取 上 述 多 传感器 数据 的 过 程 中 的 采样 时 间 间 
隔 ，S (2 ) 表示 平滑 计算 后 的 振动 信号 多 传感器 数据 ， 


大 表示 进行 平 请 处 理 


感 器 数据 的 平滑 处 理 ，; 


的 平滑 参数 。 完 成 对 振动 信号 多 传 
步 提高 数据 的 质量 。 在 上 述 


基础 上 ， 对 上 述 经 过 平滑 计算 的 多 传感器 数据 进行 融合 
处 理 0 1。 振 动 信号 多 传感器 数据 融合 的 具体 过 程 如 下 所 
示 : 


q; =m, *S(x,) 


(3) 
X= >》 9， Xi 


在 上 述 公 式 中 ，g, 表示 振动 信号 多 传感器 数据 在 不 


同 传感器 上 的 权重 值 ，m 表示 多 传感器 数据 的 加 权 值 


蕊 表示 振动 信号 多 传感器 数据 融合 的 结果 。 
1.2 基于 DBN 的 核电 机 组 振动 信号 特性 分 析 

DBN 全 称 为 深度 新 年 信念 网 络 , 是 深度 学 习 模 型 的 
一 种 ， 能 够 自 下 而 上 地 学 习 特 征 ， 对 输入 的 特征 进行 有 
效 的 特征 表示 和 分 了 0。 将 其 应 用 到 核电 机 组 振动 信号 
特性 分 析 中 ， 能 够 逐 层 提 取 核 电机 组 振动 信号 的 特征 ， 
并 将 其 进行 转化 。 基 于 DBN 的 核电 机 组 振动 信号 特性 
分 析 的 具体 过 程 如 图 1 所 示 。 


振动 信号 多 传感器 


数据 融合 结果 


提取 数据 特征 


功率 谱 质 


谐 波 


计算 特征 敏感 性 指 
数 和 模糊 粒 


输出 特性 分 析 结 果 


结束 


图 1 核电 机 组 振动 信号 特性 分 析 流 程 


如 图 1 所 示 ， 将 上 述 融 合 处 理 的 


展 动 信 号 多 传感器 


数据 作为 基础 , 利用 DBN,， 逐 层 提取 振动 信号 的 特征 如 
偏 态 特征 、 容 度 特征 、 功 率 谱 质 心 特征 、 谐 波 特征 等 
根据 提取 的 结果 ， 分 别 计算 特征 的 敏感 性 指数 和 模糊 灶 
[12] 。 
利用 DBN 提取 的 核电 机 组 振动 信号 特征 的 具体 过 程 如 
下 所 示 : 


82 =X, 三 二 


g, = RD*e dr 3 
L 


0 


在 上 述 公 式 中 ，g| 表示 振动 信号 的 偏 态 特征 ，u 表 
示 振 动 信号 的 平均 值 ，c 表示 中 位 数 ，o 表示 振动 信号 
的 方差 ，g, 表示 提取 的 裕 度 特征 ， XX, 表示 提取 的 信 


max 


号 最 大 幅 值 ，Xi, 表示 振动 信号 的 最 小 幅 值 ， 83 表示 


min 


振动 信号 的 功率 谱 质心 特征 ， 尺 (1) 表示 振动 信号 的 自 


相关 函数 ，i 表示 虚数 单位 ，g 表示 振动 信号 的 谐 波 特 


征 , 7 表示 谐 波 电流 的 波动 值 , /表示 基 波 电流 波动 值 。 


通过 计算 振动 信号 特征 的 敏感 性 指数 和 模糊 米 ， 实 现 对 
振动 信号 的 特性 分 析 中 。 其 具体 计算 过 程 如 下 所 示 : 


在 上 述 公 式 中 ， 工 表示 计算 的 振动 信号 特征 的 敏感 


了 


性 指数 ，m 表示 振动 信号 特征 的 数量 ，M 表示 振动 信 
号 特征 的 维度 ，B, 表 示 振 动 信号 的 特征 频率 值 ，, 表 


示 振 动 信号 的 相对 粹 值 ，g, 表示 振动 信 


1.3 构建 核电 机 组 振动 故障 检测 模型 

对 核电 机 组 振动 信号 进行 分 类 ， 根 据 分 类 结果 ， 实 
现 对 核电 机 组 振动 故障 的 检测 号 。 其 构建 的 核电 机 组 振 
动 故障 检测 模型 的 具体 计算 公式 如 下 所 示 : 


号 的 模糊 炳 值 。 


在 上 述 公 式 中 ， 太 ,表示 构建 的 核电 机 组 振动 故障 


检测 模型 ，o 表示 核电 机 组 振动 信号 的 序列 方差 ， 


万 (x ) 表示 核 电机 组 振动 信号 的 分 类 函数 。 完成 对 核电 


机 组 振动 信号 故障 检测 模型 的 设计 ， 将 该 检测 模型 作为 
基础 ， 对 其 进行 求解 ， 根 据 求解 结果 ， 对 其 故障 进行 检 
测 。 其 求解 过 程 如 图 2 所 示 。 


import numpy as np 

import pandas as pd 

from sklearn.model_selection import train_test_split 
from sklearn.preprocessing import StandardScaler 
from sklearn.svm import SVC 

data = pd.read_csv(Cnuclear_plant_vibration_data.csv) 
features = data[['RMS', 'Peak'", "Kurtosis']] 

labels = data['Fault'] 

X_train, X_test, y_train, y_test = train_test_split(features, labels, test_size=0.2, 
random_state=42) 

scaler = Standardscaler() 

X train = scaler.fit_ transform(X_train) 

X_test = scaler.transform(X_test) 

model = SVC(kernel='linear, C=1.0) 
model.fit(X_train, y_train) 

y_pred = model.predict(X_test) 

accuracy = np.mean(y_pred == y_test) 

print( “Accuracy: accuracy) 


图 2 核电 机 组 振动 信号 故障 检测 模型 求解 的 部 分 代码 

如 图 2 所 示 ， 设 定 阀 值 ， 判 断 核电 机 组 振动 信号 是 
否 存 在 故障 ， 并 输出 相应 的 检测 结果 [1。 
2 实验 测试 


为 提高 实验 结果 的 可 靠 性 ， 将 本 文 设计 的 方法 与 传 
统 的 方法 进行 对 比 。 其 中 ， 本 文 设计 的 基于 DBN 和 多 
传感器 数据 分 解 的 核电 机 组 振动 故障 检测 方法 为 方法 
1， 基 于 郊 狼 算法 的 检测 方法 为 方法 2， 基 于 深度 学 习 的 
检测 方法 为 方法 3。 

以 某 电力 系统 为 例 ， 以 该 电力 系统 的 核电 机 组 为 实 
验 对 象 ， 以 NVIDIA RTX2060 6G 为 显卡 ，512GB SSD 
为 硬盘 ， 在 Pycharm 为 开发 平台 ， 由 此 构建 相应 的 实验 
环境 。 其 具体 的 实验 环境 如 图 3 所 示 。 


图 3 实验 环境 
如 图 3 所 示 , 获取 某 电力 系统 核电 机 组 的 振动 信号 ， 


并 将 获取 的 振动 汇总 到 高 性 能 计算 机 中 ， 利 用 计算 机 对 

其 进行 融合 处 理 。 获 取 的 茶 电 力 系统 核电 机 组 振动 信号 
多 传感器 数据 如 图 4 所 示 。 
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图 4 某 电力 系统 核电 机 组 振动 信号 多 传感器 数据 

如 图 4 所 示 ， 即 便 是 相同 的 振动 信号 ， 其 振动 信号 
的 幅 值 也 会 存在 一 定 的 差异 性 。 

以 方法 的 AUC-ROC 曲线 为 评价 指标 ， 对比 上 述 三 
种 方法 的 性 能 。 实 验 中 ， 利 用 三 种 方法 对 某 电力 系统 核 
电机 组 的 振动 信号 进行 故障 检测 ， 统 计 其 检测 结果 的 
AUC-ROC 曲线 面积 。 其 具体 统计 结果 如 图 5 所 示 。 
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图 5 三 种 方法 的 检测 结果 

如 图 5 所 示 ， 方 法 1 的 AUC-ROC 曲线 面积 最 大 ， 
数值 更 接近 于 1， 方 法 2 和 方法 3 的 AUC-ROC 曲线 面 
积 较 小 ， 数 值 更 接近 于 0.5。 

为 进一步 验证 上 述 三 种 方法 的 性 能 ， 以 方法 的 漏 检 
率 为 评价 指标 ， 进 行 实验 测试 。 利 用 上 述 三 种 方法 对 多 
组 核电 机 组 振动 信号 进行 检测 ， 统 计 其 检测 结果 的 漏 检 
率 。 其 结果 如 表 1 所 示 。 

表 1 三 种 方法 的 检测 结果 


实验 次 数 三 种 方法 的 漏 检 率 /% 
方法 1 方法 2 方法 3 
1 0.98 10.25 25.64 
2 0.78 11.36 22.69 
3 0.96 12.54 23.54 
4 0.85 13.25 21.59 
5 0.96 13.54 23.69 
6 0.93 12.56 23.54 
7 0.85 14.58 21.56 
8 0.88 12.69 25.98 
9 0.65 13.59 26.37 
10 0.62 17.54 29.12 


如 表 1 所 示 ， 方法 1 的 漏 检 率 较 低 ， 检 测 效果 较为 
稳定 ， 漏 检 率 均 在 1% 以 下 。 方 法 2 和 方法 3 的 漏 检 率 


较 高 ， 且 数值 波动 性 较 大 。 
3 结束 语 


基于 DBN 和 多 传感器 数据 分 解 的 核电 机 组 振动 故 
障 检测 方法 ， 通 过 深度 学 习 框架 和 多 传感器 数据 融合 技 
术 ，AUC-ROC 曲线 面积 更 接近 于 1， 漏 检 率 更 低 ， 检 
测 效果 更 好 ， 有 效 提 高 了 故障 检测 的 准确 性 和 可 靠 性 ， 
能 够 从 复杂 的 多 传感器 数据 中 提取 出 与 故障 相关 的 特 
征 ， 可 以 人 全面、 准确 地 反映 机 组 的 振动 状态 ， 弥 补 了 单 
一 传感器 数据 采集 的 不 足 ， 并 利用 深度 信念 网 络 进行 分 
类 和 预测 ， 从 而 实现 快速 准确 的 故障 检测 。 
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